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Abstract 

A quasiparticle description of the thermodynamics of deconfined matter, re- 
producing both the perturbative hmit and nonperturbative lattice QCD data 
at finite temperature, is generalized to finite chemical potential. By a flow 
equation resulting from Maxwell's relation, the equation of state is extended 
from zero to non-zero quark densities. The impact of the massive strange fla- 
vor is considered and implications for cold, charge-neutral deconfined matter 
in /5-equilibrium in compact stars are given. 
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I. INTRODUCTION 

The equation of state (EoS) represents an important interrelation of state variables de- 
scribing matter in local thermodynamical equilibrium. All microscopic characteristics are 
integrated out and only the macroscopic response to changes of the state variables as, e. g., 
the temperature T and the chemical potential are retained. Via the Gibbs equation, 
e = Ts + fin — p, the energy density e is related to the entropy and particle densities. 
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s = dp/dT and n = dp/d^, respectively, and the pressure p. The thermodynamical po- 
tential p{T,fi) thus provides all information needed to evaluate, e.g., sequences of stellar 
equilibrium configurations by means of the Tolman-Oppenheimer-Volkov equation, or the 
evolution of the universe via Friedman's equations, or the dynamics of heavy-ion collisions 
within the framework of relativistic Euler equations. In the examples mentioned not only 
excited hadron matter is of relevance, rather, at sufficiently high density or temperature, a 
plasma state of deconfined quarks and gluons is the central issue. 

Within the quantum field theoretical standard model, the interaction of quarks and glu- 
ons is described by Quantum Chromodynamics (QCD). The challenge, therefore, consists 
in the derivation of the EoS of deconfined matter directly from QCD. In first attempts |l| 
the EoS of cold quark matter was derived perturbatively up to order 0{g'^) in the coupling 
g. Finite temperatures have been considered, where the perturbative expansion is nowadays 
extended up to the order 0{g^) 0. However, in the physically relevant region, due to the 
large coupling, perturbative methods seem basically to fail and, consequently, nonperturba- 
tive evaluations are needed. Present lattice QCD calculations can accomplish this task, and 
indeed the EoS of the pure gluon plasma [0,^ and of systems containing two or four light 
dynamical quark fiavors are known reliably at finite temperature, and simulations with 
physical quark masses will possibly be available in the future. 

As a matter of fact, the current QCD lattice calculations are yet restricted to zero 
chemical potential fi. Attempts to extend the lattice calculations systematically to non-zero 
chemical potential are under way |^, but focus currently more on conceptional questions 
than on precise predictions as in the case of /x = 0. However, a detailed understanding of 
the EoS at non-vanishing net quark density is required for a variety of physical questions 
as, e.g., the structure of quark cores in massive neutron stars, the baryon contrast prior 
to cosmic confinement, or the evolution of the baryon charge in the mid-rapidity region of 
central heavy-ion collisions. 

Being aware of the urgent need for this EoS, we here suggest an approach, based on a 
quasiparticle description of quarks and gluons, to map available lattice data from /i = 

2 



to finite values of /i and to small temperatures. Confidence in the model is gained by its 
quantitative agreement with the latest lattice data, which are analyzed in Sec. II, and by 
the smooth interpolation to the asymptotic limit. Discussed as well is the connection to the 
thermodynamic hard-thermal- loop resummation proposed recently p|-pl]], which provides 



a more rigorous motivation of our phenomenological approach. In Sec. Ill, we formulate 
the extension of the model to finite chemical potential, which does not require any further 
assumption. Although our approach is formulated in terms of gluon and quark quasiparticles 
and, hence, designed for 'conventional' deconfined matter, there are formal indications of 
the existence of other QCD phases. In Sec. IV the quark-gluon plasma is considered with 
the massive strangeness degree of freedom appropriately taken into account. Despite of two 
free model parameters (which can be fixed when lattice data for the physical (2+1) fiavor 
system will be available), our approach makes some definite statements about the EoS of the 
quark-gluon plasma. This is shown in a detailed study of the thermodynamics of /^-stable 
deconfined matter, which allows us to draw conclusions about the bulk properties of pure 
quark stars. Finally, a summary and an outlook are given in Sec. V. 

II. THE QUASIPARTICLE MODEL 

The concept of describing an interacting system in terms of a system of quasiparticles 
which, as appropriate ground states, refiect relevant aspects of the interaction by effective 
properties, is successfully used in many fields of physics. For deconfined QCD plasmas, the 
development of the quasiparticle idea, starting with [p!2|-[T^ , ran parallel to the improvement 
of the lattice QCD calculations. These nonperturbative numerical simulations provide the 
testing ground for the quasiparticle concept which is a phenomenological approach. At the 



same time, by identifying the quasiparticles with the asymptotically relevant excitations [15 



these models yield a physically intuitive picture of the strongly interacting system. Referring 



to |T6[ for a recent survey of quasiparticle descriptions of finite-temperature lattice data, we 



focus in this work on the generalization of the approach to non-zero chemical potential. 
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A. Formulation of the model 



We consider an SU{Nc) plasma of gluons, with Nc = 3 for QCD, and Nf quark flavors 
in thermodynamic equilibrium. Within the approach outlined below, the interacting plasma 
is described as a system of massive quasiparticles, a picture arising asymptotically from the 
in-medium properties of the constituents of the plasma. For thermal momenta A; ~ T, the 
relevant modes, transversal gluons and quark-particle excitations, propagate predominantly 
on-shell with dispersion relations ujf{k) ~ -|- k"^ and 



(1) 



while the plasmon (longitudinal gluon) and plasmino (quark-hole) states are essentially un- 
populated [|17|. Neglecting sub- leading effects, the medium contributions 11* are the leading 
order on-shell selfenergies. Depending on the coupling G^, the temperature and chemical 
potential as well as the rest mass moi (the latter vanishing for gluons), the 11* are given by 
the asymptotic values of the gauge independent hard-thermal/density-loop selfenergies |p!7| . 
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We emphasize that the contributions IT* are generated dynamically by the interaction within 
the medium, so the quantities in (|I]) can be considered only as effective masses, with no 
ambiguities arising for the gauge boson quasiparticles. 

Generalizing the approach of Ref. to a finite chemical potential fi controlling a con- 
served particle number, the pressure of the system can be decomposed into the contributions 
of the quasiparticles and their mean field interaction B, 
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where pi = ±diT f d?k/{2TiY ln(l ± exp{ — (cjj — jj,i)/T}) is the pressure of an ideal gas of 
bosons or fermions with degeneracy rfj and the state dependent effective masses (|l],|D. By the 



general stationarity of the thermodynamic potential Q = —pV under functional variation of 



the selfenergies |T9|, which in the present approach simplifies to dp/dH* = 0, i? is related 
to the quasiparticle masses, 

OB ^ dp,{T,fif,m]) 

dU* dm] ■ ^ ^ 

We remark that the stationarity of the quasiparticle pressure leads, by the Feynman- 



Hellmann relation, to a chiral condensate {ipip) ~ ^oq, so the restoration of chiral symmetry 
of massless quark flavors in the deconfined phase is inherent in the approach. Furthermore, 
the stationarity implies that the entropy and the particle densities are given by the sum of 
the quasiparticle contributions. 



_ dpi{T, jj,i;m; 

Si 



dT 



, -"'^^ ^ . (5) 
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This quasiparticle approach reproduces the leading-order perturbative results for arbitrary 
values of Nc and Nf. Moreover, however, it represents a thermodynamically consistent 
effective resummation of the leading-order thermal contributions . As a nonperturbative 
approach, hence, it can be expected to be more appropriate in the large-coupling regime 
than perturbative predictions which, as truncated power expansions, fail to be reliable there. 
More specifically, the effective quasiparticle description will be an appropriate framework 
as long as the spectral properties of the relevant excitations do not differ qualitatively 
from their asymptotic form. Since propagating constituents in a plasma are, on general 
grounds, expected to possess a spectral representation with the strength predominately 
accumulated in a quasiparticle peak, our approach may be anticipated to be reasonable 
even close to the confinement transition. In a way, this expectation has some support from 
the hot 0^ theory where the selfconsistent resummation of superdaisy graphs leads to a 
similar 'strict-quasiparticle' picture |21|] as introduced above for QCD. By resumming the 



propagator beyond the leading-loop order, then, the quasiparticle properties were shown 
to be pronounced even for larger values of the coupling In a more direct way, the 

quasiparticle approach can be tested by comparison to the lattice data available for several 
QCD-model systems at finite temperature. 



B. Test of the model 



To obtain the EoS as a function of the temperature, G'^{T) has to be specified in Eq. (^. 
The quantity is to be considered as an effective couphng since it parameterizes all de- 
viations of the exact spectral function from our 'strict-quasiparticle' ansatz. On the other 
hand, to make contact to perturbation theory, it has to approach the running coupling in 
the asymptotic limit of large temperatures. In the form 

G\T^^^ = ^) = ^ (6) 

{llNc-2Nf)hi' " 



TJX 

the shape of the exact spectral function is encoded in the single number Tg/Tc-, whereas Tj A 
is related to the QCD scale Aqcd- While the physically intuitive ansatz turns out to be 
applicable to all systems considered in the following, other parameterizations are conceivable 
as well 



As a first example, the quasiparticle model is applied to the case of the pure SU(3) gauge 
plasma. Recently, new lattice results calculated with a renormalization-improved action have 
been published [Q, which are similar but by some 3 ■ ■ ■ 4% larger than the data 0] obtained 
with a standard action. Consequently, the parameters fitting both data sets, as shown in 
Tab. 1, differ notably only in the number dg of the gluon quasiparticle degrees of freedom. 
The quantity dg is introduced here as fit parameter to account, on one hand, for possible 
finite size effects of the lattice data and, on the other hand, for residual sub-leading effects 
otherwise neglected in the quasiparticle ansatz. As expected, for both data sets the values 
of dg are close to 2{Nl — 1) = 16. As shown in Fig. |l] for the entropy, which is the quantity 
to adjust the parameters according to Eq. (j^), the quasiparticle model reproduces the lattice 
results quantitatively. In fact, the deviations are smaller than the numerical uncertainty of 
the data. 

The next case analyzed is the system with two light quark flavors, with an estimate 
for the continuum extrapolation of the pressure given in [0]. To calculate the quasiparticle 
pressure @ for fixed values of the parameters, the interaction pressure B is required as a 



function of the temperature. It can be obtained by the relation (H), with an integration 
constant Bq = B{Tq) adjusted to the lattice data, e.g., at the smallest temperature. The 
resulting agreement of the data and the model, with the parameters given in Tab. 1, is shown 
in Fig. 1^. Note that the quark quasiparticle degeneracy is held fixed to dq = 2{Nf-i) '^f 
in the next example. 

The available lattice data |P for the system with Nf = 4 light flavors were obtained by 
the extrapolation to the chiral limit of lattice simulations of quarks with a ratio of mass to 
temperature of 0.2 and 0.4. Being compatible with each other, both data sets are reproduced 
by a single set of model parameters as shown in Fig. ^. So far, however, the data are not yet 
continuum extrapolated. Accordingly, the quasiparticle degeneracies are considerably larger 
than in the previous cases (cf. Tab. 1). Thus it is suggesting to interpret, to a large amount, 
the fit value of dg/16 to be due to finite size effects which, in the free limit, are indeed found 
to have the same order of magnitude. Concerning the function B{T) we mention that it 
turns out positive in the interval 1 < T/T^ < 2, and negative and small at larger values of 
T 0. 

We conclude this section with two remarks. As anticipated, the quasiparticle approach 
(being tested with available lattice data) is an appropriate description even close to the 
confinement transition. This provides some confidence in the applicability, in the nonper- 
turbative regime, of the phenomenological quasiparticle resummation with the effects of 
finite spectral widths parameterized in an effective coupling. The quasiparticle approach is 
also supported by the more formal development of the thermodynamical resummation of 
hard-thermal loops, proposed in to improve the bare perturbation theory. As shown in 
T0| , the entropy of the plasma is dominated by the asymptotic behavior of those excitations 



which we start with as the relevant quasiparticles. Moreover, a direct hard-thermal-loop 
resummation of the pressure within the Luttinger-Ward formalism results in an expression 
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similar to Eq. (|^), with B corresponding in part to the so-called $ functional] [^. While 
this resummation, being more rigorous from a theoretical point of view and containing no fit 
parameters, reproduces the lattice data already qualitatively, we focus here on a phenomeno- 
logical description, adjustable quantitatively to data. Nevertheless, the hard-thermal- loop 
resummation approaches yield further confidence to the quasiparticle model, also with regard 
to the extension to finite chemical potential as outlined subsequently. 

III. EXTRAPOLATING TO FINITE CHEMICAL POTENTIAL 

Encouraged by the successful quasiparticle description of the = lattice data, the 
model is now applied to finite chemical potential. This extension requires no further as- 
sumptions and allows in a straightforward way to map the EoS at finite temperature and 
fi = into the fiT plane. This continuous mapping relying on quark and gluon quasi- 
particles, however, cannot provide information about other possible phases with a different 
(quasiparticle) structure, so the following consideration apply only to 'conventional' decon- 
fined matter. 

In general, the pressure is a potential of the state variables T and fi. As a direct conse- 
quence thereof, the Maxwell relation implies for the quasiparticle model 

E 

i 

which, formally, is the integrability condition for the function B defined by Eq. (^. With 
n* depending on G^, Eq. (0) represents a fiow equation for the effective coupling, following 
directly from principles of thermodynamics. This flow equation is a quasilinear partial 
differential equation of the form 



^ Following these parallels further, we note that the hard-thermal-loop resummed entropy is 
independent of <I>, as the quasiparticle entropy Eq. (|5|) is independent of B. This feature, related to 
the underlying one-loop structure of the selfenergies in both cases, makes the entropy an interesting 
quantity, whereas the pressure contains the full thermodynamical information. 
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drij dU* dsj dll* 
dmf dT dmf djj, 



(7) 



with the coefficients a^,^ and b depending on T, /i and G^. It is instructive to consider ffist 
the asymptotic hmit, where G^ ^ <^ 1 and Eq. reduces to the homogeneous form 



with c = (4A'"c + 5Nf) /{9Nf). This equation has solutions = const along the charac- 
teristics given by cT^ + 2T^(yu/7r)^ + (/u/vr)"^ = const. As a result of this elliptic flow, the 
renormalization scale Tc/A of g^{T,0) determines the scale TcTtc^^^/X of the /x dependent 
running coupling g^{0,fi). 

The flow of the effective coupling is elliptic also in the nonperturbative regime, thus 
mapping the /i = axis, where G^ can be determined from lattice data, into the /iT 
plane. As a representative example, the nonperturbative flow of the coupling of the Nf = 4 
system is shown in Fig. ^ As to be expected, the characteristics emanating at /i = 
in the deconfined phase do not attach a certain region, at small values of T and fi, of 
the phase diagram which is obviously related to the confinement phase. It is noted that 
at /i ^ 3Tc and small temperatures, indicated by intersecting characteristics, the solution 
of the flow equation is not unique. This ostensible ambiguity is, however, of no physical 
consequence, i.e., the approach is intrinsically consistent, since already outside that region 
the quasiparticle pressure (calculated along the characteristics) would become negative. This 
region of absolute thermodynamical instability, on the other hand, may be considered as an 
indication of a transition to a different phase, in fact already at some positive pressure. 
Whether this is the confined phase or, maybe already at such small chemical potential, a 
color-superconducting state (cf. p4[ and further references therein) remains, of course, a 
speculation within the present phenomenological approach. In any case, the pressure (see 
Fig. ^) increases much slower with /i than with T, which makes the transition to another 
deconfined phase likely in that region. This behavior is mainly attributed to the fact that 
the effective coupling (see Fig. |^), as a function of increasing chemical potential, approaches 
the asymptotic limit more slowly than in the direction of increasing temperature. 



To conclude this example, we remark that the proposed extension of the lattice EoS to 
finite chemical potential, though mo del- dependent in using a quasiparticle picture to derive 
the fiow equation (^), is based on ab initio calculable data. Since the model works perfectly 
at yU = and, on the other hand, the convergence problem of perturbation theory and the 
concept to overcome it by resummation are not specific to /i = 0, we expect our extension of 
the data as reliable, at least semi-quantitatively. Similar fiow equations can also be derived 
in more rigorous approaches, and on general grounds they are expected to be also of elliptic 
type. In this sense, the EoS at ;U = 'knows' about 7^ 0. Nevertheless, as already 
mentioned, such continuous extensions do not access other possible Q CD-phases with a 
different ground state; the actual EoS of strongly interacting matter is rather determined 
by thermodynamical stability between 'competing' phases. 



IV. APPLICATION OF THE MODEL 

A. Strangeness in matter 

Although lattice simulations of the physically interesting case of two light and one 
medium- mass fiavor are in progress reliable data of the thermodynamics of QCD with 
physical quark current masses, i. e. with mps ~ 150 MeV for the strange fiavor, are still 
lacking to fix the quasiparticle parameters at /i = 0. Nevertheless, even when freely varying 
the parameters over large ranges, our model allows some predictions since the parameter 
dependence of the quasiparticle pressure p'^^ is weak due to the stationarity property Eq. (||). 
The parameters are restricted to match p'^^ to the hadronic pressure p^"''^ at the confinement 
temperature which we assume to be Tc = 150 MeV. The uncertainty of p^"'^(Tc, 0) = 3.1 ■ 10^ 
MeV^, as estimated by a hadron resonance gas model, can be absorbed into the integra- 
tion constant Bq = B{Tc,0). This parameter, being related phenomenologically to the bag 
constant, is varied in the range 120 MeV < -B^^ < 180 MeV. (Larger values of Bq yield an 
increased ratio e/T^ close to Tc, opposed to what is expected from available lattice data.) 
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For the second independent parameter in p*^^, the range 3 < A < 11 is considered reasonably 
large, as suggested by the parameters of the analyzed lattice data in Tab. 1. The EoS for a 



representative choice of the parameters is provided in tabular form in |^6[| for the plasma with 
vanishing net-strangeness, i.e. with the quark chemical potentials /is = and fiu = fJ'd = A*- 
With the effective coupling decreasing only logarithmically from a = G"^ j (47r) ~ 1 at small 
values of the temperature and the chemical potential, the medium contributions (H) to the 
quasiparticle masses are large, of the order of T or /i. Hence, the quasiparticles are consid- 
erably under-saturated in phase space compared to the asymptotic limit, while the current 
mass of strange quarks turns out a less important suppression factor. 

B. Deconfined /3-stable strange matter 

With the aim to study implications for quark matter stars, we consider in the following 
a plasma of gluons, quarks and leptons in /3-equilibrium maintained by the reactions ci, s 
u + / + I//, which imply the relations /id = /^s = A^m + /i« = among the chemical potentials. 
The lepton chemical potential /i/ is determined as a function of T and /i by the requirement 
of electrical charge neutrality. Without noticeable change of the result, the lepton pressure 
can be estimated by either the pressure of a free electron gas or, doubling the degrees 
of freedom, a gas of free electrons and muons. The pressure p = + ]} and the resulting 
energy density at /i = are shown in the left panels of Fig. ^ for several values of the 
parameters A and i?o- Already at temperatures slightly above Tc, the scaled energy density 
reaches a saturation-like behavior at about 90% of the asymptotic limit. This qualitative 
feature, known from the lattice simulations 0-0], is to a large extent insensitive to the 
specific choice of the free parameters, while Bq has a distinct impact on the latent heat. The 
resulting EoS at vanishing temperature is displayed in the right panels of Fig. |^. In this 
regime, the asymptotic values are approached more slowly due to the less rapid decrease of 
the effective coupling with increasing values of /i, similar to the N f = 4 plasma. For various 
sets of parameters, the EoS is available in tabular form PB . 
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C. Pure quark stars 



By the Tolman-Oppenheimer-Volkov (TOV) equation (cf. Eq. (2.212) of Ref. [^), se- 
quences of hydrostatic equihbrium configurations of stars can be calculated, given the rela- 
tion e{p) of the star matter. Since the bulk properties of the star are less dominated by the 
outermost shells, we consider the case of pure quark stars as a reasonable approximation 
for compact stars with a large quark matter core. For energy densities up to several times 
the nuclear saturation density and at temperatures less then some 10 MeV, the dependence 
e{p) of /5-stable quark matter, as estimated numerically by our quasiparticle model, can 
be parameterized by e = AB + dtp. While the naive bag model EoS had a = 3, our EoS, 
with the considered choices of the model parameters, yields 3.1 < a < 4.5, indicating the 
nontrivial nature of the interaction. For parameters A > 5, values of B^^^ > 200 MeV are 
found, while only the extreme choice of A = 3 allows B^^^ as small as 180 MeV. However, 
the value of B, i. e., the energy density at small pressure, has a strong impact on the star's 
mass and radius, obtained by integrating the TOV equation. 

In Fig. D, the mass as function the radius of pure quark stars is displayed for various values 
of the parameter B and a, while in Fig. ^the maximum mass and the corresponding radius 
are shown as function of B^^^. For B^^^ > 200 MeV, as turned out to be the most likely range 
according to our analysis, the quark stars have masses less than one solar mass Mq. Compact 
quark stars with a larger mass can only exist for values of 5^/^ < 180. ..200 MeV, depending 
only marginally on a. We hence conclude that compact stars with masses M ~ Mq and 
radii R ~ 10 km are unlikely to be composed mainly of deconfined matter in /3-equilibrium, 
irrespective of the uncertainty of the model parameters of our approach. 

V. SUMMARY AND OUTLOOK 

In summary we have generalized a thermodynamic quasiparticle description of deconfined 
matter to finite chemical potential fi not yet accessible by present lattice calculations. Of 



12 



central importance to the model is the effective coupling G'^{T,ij) which can be obtained 
at = from available lattice data, proving at the same time the applicability of the 
quasiparticle picture even close to the confinement transition. At 7^ 0, is determined 
by a flow equation following directly from the Maxwell relation of thermodynamics. As a 
result of the elliptic flow, the basic features of the EoS at // = 0, namely the nonperturbative 
behavior near confinement and the asymptotic behavior, are mapped into the 11 T plane as 
exemplified by the Nf — A flavor system studied on the lattice at /j, — 0. 

An important consequence of the quasiparticle approach is the relation of the 'critical' 
values of temperature and chemical potential. For deconfinement matter with physical quark 
masses, this fact leads to the implication that pure quark stars are less massive and, hence, 
more distinct in the bulk properties to neutron stars than estimated by other approaches 
with an ad hoc choice of small values of the bag constant. 

Further applications of the quasiparticle model, in connection with realistic thermody- 
namical descriptions of other QCD phases, is desirable. In particular, in the cold and dense 
regime where a color-superconducting state is expected to be important, the quasiparticle 
EoS can provide one part of the information about the phase boundary. With regard to 
the on-going heavy-ion program at the relativistic heavy-ion collider RHIC at Brookhaven 
National Laboratory, we finally would like to point out the applicability of our quasiparticle 
EoS as an input to hydro-relativistic simulations. For a realistic description of the (near-to) 
equilibrium states of a heavy-ion collision, detailed knowledge about the EoS at non-zero 
baryon densities is required. Since the results of the quasiparticle model difi^er considerably 
from that of simpler approaches as the bag model (e.g., in the reduced latent heat and the 
saturation-like behavior of the energy density at about 90% of the free limit), dynamical 
effects of a realistic EoS may even lead to observable signals of the experimental formation 
of the quark-gluon plasma. 
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TABLES 





A 


Ts/T, 


dg 


SU(3) i 


5.02 


-0.75 


16.9 


SU(3) i 


4.83 


-0.72 


17.5 


Nf = 2§ 


10.2 


-1.00 


17.0 


Nf = 4§ 


6.59 


-0.80 


20.6 



TABLE I. Quasiparticle-model parameters adjusted to the lattice data. For the systems con- 
taining quarks, the quark degeneracy is fixed to dq = 2(n^^~i) ^9' 
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FIG. 1. The comparison of the quasiparticle entropy density with the lattice results, calculated 
with a standard Wilson action |^,^ and a renormalization-improved action , of the pure gluon 
plasma. 
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FIG. 2. The pressure in our model compared to the estimate of the continuum extrapolation 
of the lattice data for the light 2-flavor case. 
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FIG. 3. The comparison of the energy density in our model with the chiral extrapolation of the 
lattice data IgI of the 4-flavor case. 
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FIG. 4. The characteristics of the coupling flow equation ^ for the QCD plasma with N ^ = 4 
light flavors for which G'^{T,ii = 0) is obtained from the lattice data j^. At leading order, the 
characteristics are curves of constant coupling strength. The region of intersecting characteristics 
is of no physical relevance since the pressure becomes negative in the region below the dash-dotted 
line (see text). 
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FIG. 5. The pressure of the Nf = 4 plasma in the chiral Umit as a function of /x and T, scaled 
by free limit. 
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FIG. 7. Total pressure and energy density (lower and upper bands, respectively) of the 
charge-neutral quark-gluon plasma in /^-equilibrium, scaled by the values of the free limit. The 
panels on the left (right) show the EoS at ^ = (T = 0), for values of the model parameter 
3 < A < 11 (lower and upper line, respectively, of the shaded area), and B^'^ = 120 MeV (top) and 
bI'^ = 180 MeV (bottom). Non -unique values of the energy only occur in the unphysical region, 
where p < 0. 
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FIG. 8. The dependence of the mass of pure quark stars on the radius for several values of the 
parameters a and B^^^. 
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FIG. 9. The dependence of the maximum quark star mass (upper band and right scale) and 
the corresponding radius (lower band and left scale) on the parameter B^^^. The upper (lower) 
boundary of the hatched bands are for a = 3 (4.5). 
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